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Abstract 

The Large Hadron Collider forward (LHCf) experiment is designed to use the LHC to verify the hadronic-interaction 
models used in cosmic-ray physics. Forward baryon production is one of the crucial points to understand the devel¬ 
opment of cosmic-ray showers. We report the neutron-energy spectra for LHC \[s -1 TeV proton-proton collisions 
with the pseudo-rapidity t) ranging from 8.81 to 8.99, from 8.99 to 9.22, and from 10.76 to infinity. The measured 
energy spectra obtained from the two independent calorimeters of Arml and Arm2 show the same characteristic 
feature before unfolding the difference in the detector responses. We unfolded the measured spectra by using the 
multidimensional unfolding method based on Bayesian theory, and the unfolded spectra were compared with cur¬ 
rent hadronic-interaction models. The QGSJET 11-03 model predicts a high neutron production rate at the highest 
pseudo-rapidity range similar to our results and the DPMJET 3.04 model describes our results well at the lower 
pseudo-rapidity ranges. However no model perfectly explains the experimental results in the whole pseudo-rapidity 
range. The experimental data indicate the most abundant neutron production rate relative to the photon production, 
which does not agree with predictions of the models. 
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1. Introduction 

The forward particle production process induced by 
collisions of high-energy particles is a poorly under¬ 
stood phenomenon in high-energy physics. Though it 
is important to understand the development of cosmic- 
ray showers in the atmosphere, the validity of hadronic- 
interaction models have not been sufficiently verified at 
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energies for ultra-high-energy cosmic rays (UHECRs, 
> 10 18 eV) because of lack of experimental data in this 
energy range. This lack of data results into a large un¬ 
certainty in the interpretation of the energy and chem¬ 
ical composition of UHECRs. Forward baryons play a 
very important role in the development of cosmic-ray 
showers. If forward baryons carry more collision en¬ 
ergy, cosmic-ray showers develop much deeper in the 
atmosphere, and vice versa. However, in the energy 
range of UHECRs, the predictions by current models 
differ significantly among themselves. 
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The excess of muons at ground level is reported as 
one of the problems in the cosmic-ray shower observa¬ 
tions. The number of muons observed by the surface 
detector array of the Pierre Auger Observatory (PAO) 
fill is higher than the one expected based on the en¬ 
ergy determined by the fluorescence detectors even if 
a heavy primary mass is assumed 0. It is suggested 
that the number of (anti) baryons generated in the for¬ 
ward region is strongly related to the number of muons 
observed by PAO at the ground 0. Therefore, baryon 
production at the very forward region is quite important 
to understand cosmic-ray showers. 

In this paper, we report the results of analyzing the 
data of the Large Hadron Collider forward (LHCf) 
experiment for forward neutron spectra. Forward 
baryon spectra at the laboratory equivalent energy of 
2.5xl0 16 eV ( V^=7TeV) will be a crucial input to im¬ 
prove the hadronic-interaction models used in the air 
shower analyses. 

2. LHCf experiment 

The LHCf experiment is designed to use the LHC to 
verify the hadronic-interaction models used in cosmic- 
ray experiments 0 i. Two independent detectors 
named Arml and Arm2 were installed in the detector 
installation slots of the TANs located 140 m away from 
the interaction point 1 (IP1). Because charged parti¬ 
cles are swept away by the D1 bending magnets, LHCf 
can measure only neutral particles in the very forward 
region of the LHC (pseudo rapidity I 77 I > 8.4). Both 
detectors have two different sampling calorimeters with 
44 radiation lengths (1.6 hadron-interaction lengths) of 
tungsten plates and 16 layers of sampling scintillators 
0]. Four layers of the position sensors (SciFi in Arml 
and silicon micro-strip sensor in Arm2) can measure the 
hit position transverse to the beam direction. The trans¬ 
verse dimensions of the calorimeters are 20 mm x 20 
mm and 40 mm x 40 mm in Arm 1, and 25 mm x 25 mm 
and 32 mm x 32 mm in Arm2. The details of the detec¬ 
tor performance during the 2009-2010 proton-proton 
collisions are reported in 0 . 

The performance of the LHCf detectors for hadron 
measurements was studied by Monte Carlo (MC) simu¬ 
lations and confirmed by using 350 GeV proton beams 
at CERN-SPS 0. Depending on the incident-neutron 
energy, the energy resolution and position resolution are 
about 40% and 0.1-1.3 mm, respectively. The detection 
efficiency for neutrons was estimated to be 70%-80% 
for neutrons above 500 GeV. 

In this paper we assume hadronic showers are pro¬ 
duced by neutrons. According to the EPOS 1.99 gener¬ 


ator 0, about 10% of other hadrons, i.e.. As and K°s, 
are also included in the data and this fraction and the 
energy dependence are model dependent. 


3. Analysis 


3.1. Data used in analysis 

The data used in this analysis were obtained on May 
15, 2010 from proton-proton collisions at y[s — 7 TeV 
(LHC Fill # 1104). The typical luminosity correspond¬ 
ing to this fill derived from the counting rate of the 
LHCf front counters 0 was ( 6 .3-6.5) xl0 28 cnT 2 s _1 . 
The data set was the same as one used in the previously 
published photon analysis results, and additional de¬ 
tails can be found in [ id]. The trigger for LHCf events 
was generated when signals from any three successive 
scintillation layers in any calorimeter exceeded a pre¬ 
defined threshold (typically 130 minimum ionizing par¬ 
ticles (MIPs)). Data acquisition (DAQ) was performed 
under 85.7% (Arml) and 67.0% (Arm2) average live- 
times. 

Taking the DAQ live-times into account, the inte¬ 
grated luminosities of the data set were 0.68 11 b for 
Arml and 0.53 nb 1 for Arm2, each with +6.1% uncer¬ 
tainty. The numbers of inelastic collisions were about 
48M and 38M collisions for Arml and Arm2, respec¬ 
tively. 

MC predictions were conducted with the generators 
DPMJET 3.04 HI, EPOS 1.99 0, PYTHIA 8.145 
[H, QGSJET 11-03 fH, and SYBILL 2.1 Q and 
compared with the experimental results. In the MC sim¬ 
ulations, the COSMOS (v7.49) and EPICS (v8.81) fl]] 
libraries that are used in air-shower and detector simu¬ 
lations were used to simulate the flight of particles from 
the IP1 to the detectors and the response of the detec¬ 
tors. About 10M inelastic collisions were simulated for 
each model. 


3.2. Event reconstruction 

Initially, the offline-event selection was applied when 
energy depositions equivalent to more than 200 MIPs 
were recorded for three successive layers in addition to 
the experimental trigger. The positions at which parti¬ 
cles hit the detector were determined by using the po¬ 
sition sensors. Because reconstruction of the events is 
difficult at the edge of the calorimeters due to large fluc¬ 
tuations in the energy deposition, events within 2 mm 
from the edge were discarded from the analysis. The 
lateral shower leakage caused by the limited lateral size 
of the detectors deteriorates the energy resolution. This 
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Figure 1: The L 2 D parameter distribution for the experimental data 
and the MC simulations from the template MC. The closed circles 
represent the Arml experimental results, whereas the red and blue 
histograms correspond to photon and neutron predictions. The open 
circles is the scaled results of the MC simulation obtained by method 
A. 


position-dependent leakage effect was corrected by us¬ 
ing the transverse-hit position measured by the position 
sensors. 

Particle identification (PID) between neutrons and 
photons was based on the difference in the longitudinal 
shape of the shower development. Two simple param¬ 
eters called Z. 20 % and Lgo% were introduced to charac¬ 
terize the shower shape. These parameters were defined 
as the depths containing 20% and 90%, respectively, of 
the total energy deposited within the layers. Consider¬ 
ing the correlation between Z/>o% and Lgo% an optimized 
parameterL 2 D was defined as L 2 D = Lgo% - l/4x £ 20 % 
to improve the selection efficiency and purity than previ¬ 
ous analyses. FigureQ]shows the L 2 D distributions. Two 
distinct peaks are identified in the observed L 2 D distri¬ 
bution indicated by the closed circles. The histograms 
correspond to the MC prediction of pure photons (red) 
and pure neutrons (blue) and they are called ‘templates’ 
hereafter. The templates were produced by accumulat¬ 
ing the MC simulation as a mixture of five models, DP- 
MJET 3.04, EPOS 1.99, PYTHIA 8.145, QGSJET IT- 
03, and SYBILL 2.1, with same statistics for each. 

To obtain neutron spectra, only events with the pa¬ 
rameter L 2 D exceeding a certain threshold were identi¬ 
fied as neutron-like events. The effects of the neutron 
selection efficiency and the photon contamination were 
corrected by using the efficiency e and purity P (i.e.. 


multiplying by P/e ) determined by MC simulations 
and the template fitting method fiMEi, respectively. 
To estimate the photon contamination, the templates for 
photons and hadrons were independently scaled to re¬ 
produce the experimental results (called method A). The 
open circles in Figure [I] represent the fitting result with 
method A. Neutron selection criteria, as indicated in 
FigureQ] were chosen to maximize ex P. To cope with 
energy dependence, events were separated into eight 
different energy ranges according to reconstructed en¬ 
ergy. 

After correcting the PID efficiency and purity, we ob¬ 
tain the production rate of neutrons as a function of ob¬ 
tained energy that was determined from the total de¬ 
posited energy in the calorimeter and from the iden¬ 
tity of the particle. Details of the event reconstruction 
for neutrons are summarized in 0|. About 0.3 million 
neutron-like events passed the PID selection for each 
arm. 

To combine the results of Arml and Arm2, we se¬ 
lected events that occurred within the common rapid¬ 
ity regions. Events within 6 mm (77 > 10.76) from the 
beam center were selected for the small towers. The 
large towers of Arml and Arm2 were divided into two 
regions. The inner region ”A” was defined by a radius 
of 28-35 mm (8.99 < 77 < 9.22) from the beam center, 
whereas the outer region ”B” was defined by a radius 
of 35-42 mm (8.81 < 77 < 8.99). For the analysis, we 
used azimuthal-angle intervals d(f> of 360° for the small 
towers and 20° for the large towers. 


3.3. Systematic uncertainties 

Energy scale 

In order to determine energy scale and estimate its 
systematic uncertainty, we followed the previous 
analyses HoHH. We observed invariant mass excesses 
of 8.1% (Arml) and 3.7% (Arm2) compared with the 
7 r° mass reconstructed in the MC simulations. These 
excesses were corrected in this analysis. Based on 
this mass excess correction and known calibration 
uncertainty, totally +5.6% (Arml) and +4.4% (Arm2) 
were assigned as systematic uncertainties with respect 
to the central value of the mass shift. In addition, 
according to the differences between the SPS beam 
test and MC simulation in the reconstructed energy 
of 350 GeV proton showers |0], +2.0% (Arml) and 
-3.8% (Arm2) errors were quadratically added to the 
respective energy scale uncertainties. 

PID 

Method A described in section 13.21 did not perfectly 
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Table 1: Relative differences in eight bin energy spectra to evaluate systematic uncertainties from PID process. Small tower, Large tower A, and 
Large tower B correspond to the rapidity ranges from 10.76 to infinity, from 8.99 to 9.22, and from 8.82 to 8.99, respectively. 


LHCf Energy 
[GeV] 

100 

-500 

500 

-1000 

1000 

-1500 

1500 

-2000 

2000 

-2500 

2500 

-3000 

3000 

-3500 

3500< 

Small tower 

-8.0% 

-3.4% 

-1.8% 

-1.1% 

0 .0% 

0 .0% 

0 .0% 

-0.7% 

Large tower A 

-12.8% 

-9.7% 

-4.2% 

-1.7% 

-0.7% 

-0.6% 

0 .1% 

-0.8% 

Large tower B 

-17.7% 

-11.5% 

-3.6% 

-1.2% 

-0.7% 

-0.6% 

0.4% 

-1.3% 


reproduce the experimental results. To estimate sys¬ 
tematic effects from these differences, we used a more 
artificial method (Method B) that allows longitudinal 
displacements and modifications in width of the distri¬ 
butions to fit the experimental results. The systematic 
uncertainty from the PID process was estimated by 
comparing the results using the Method A and Method 
B. The energy-dependent PID systematic uncertainty 
is mostly 1% above 1.5 TeV and 12% at most below 
this energy. The relative differences in the neutron 
production rate defined as 


1 


Pb/cb 


0 ) 


Pa/ca 

are summarized in Table [7] Here, e,\ (e#) and Pa (Pb) 
are the efficiency and purity determined by Method A 
(Method B), respectively. Final results are given using 
the Method A, while the differences shown in Table [I] 
are taken as a part of the systematic uncertainties. 

Multi-hit 

When two or more particles enter one of the LHCf 
calorimeters, these events are called ‘multihit’ events. 
Because discriminating between single and multihit 
events for neutron-like events was difficult because 
of large fluctuation of hadronic showers, rejection of 
multihit events causes a large systematic uncertainty. 
Because the multihit event rate was predicted by MC 
simulations to be less than a few percent, all events 
used for this analysis were treated as single-hit events. 
The systematic effects on the energy spectra from the 
multihit events were studied by using MC simulations. 
We tested the difference of the spectra with current 
method and those with ideal multihit reconstruction by 
using the MC study. The difference was less than 6% 
above 500 GeV as summarized in Tableland was taken 
into account as a part of the systematic uncertainties. 

Position resolution 

Because the resolutions of transverse hit position for 
neutrons are 0.1-1.3 mm depending on the neutron 
energy 0], some events were reconstructed as hits at 


positions far from the true hit position. Thus, particles 
hitting outside of the fiducial area may be contaminated 
and vice versa. The effect of the position resolution was 
estimated by using MC simulations. 

The difference between neutron-energy spectra 
selected by the reconstructed position and the true hit 
positions were calculated with EPOS 1.99, QGSJET 
11-03, SYBILL 2.1, DPMJET 3.04, and PYTHIA 
8.145. Because no significant dependence on model 
was found, the average of the predictions by the 
five different models was assigned to the systematic 
uncertainty. The systematic error from these effects 
was less than 8.4% at energies above 500 GeV and are 
summarized in Table [3] 

Other systematic errors 

As similar as the previous study Ei, the other system¬ 
atic errors, such as the integrated luminosity (±6.1%) 
and position of beam center (typically ±3-10%) were 
taken into account. The systematic uncertainty from 
pile-up events (0.2%) was negligibly small. 


4. Results 

4.1. Measured energy spectra 

Figure [2] shows the energy spectra of forward neu¬ 
trons measured by the LHCf Arml and LHCf Arm2 de¬ 
tectors. The energy scale correction described above 
was applied in these spectra (-8.1% for Arml and - 
3.7% for Arm2). The vertical axes were normalized to 
the number of inelastic collisions per GeV. The vertical 
bars represent the statistical (they are very small) and 
systematic uncertainties except the energy scale and lu¬ 
minosity uncertainties Q. The closed and open circles 
show the results of the Arml and Arm2 detectors, re¬ 
spectively. A small difference in the detection efficien¬ 
cies of the Arml and Arm2 detectors were not corrected 


'The uncertainty from the luminosity determination is not indi- 
cated in the figures throughout this paper. 
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Table 2: Relative differences of ten bin energy spectra to evaluate systematic error from multihit events. Small tower, Large tower A, and Large 
tower B correspond to the rapidity ranges from 10.76 to infinity, from 8.99 to 9.22, and from 8.82 to 8.99, respectively. 


LHCf Energy 
[GeV] 

100 

-500 

500 

-1000 

1000 

-1500 

1500 

-2000 

2000 

-2500 

2500 

-3000 

3000 

-3500 

3500 

-4000 

4000 

-5000 

5000 < 

Small tower 

18.0% 

5.6% 

1 .8% 

0.1 % 

0.1 % 

0.3 % 

-0.1% 

-2.5% 

-2.2% 

-3.6% 

Large tower A 

9.5 % 

5.1% 

-0.4% 

0.4% 

0 .2% 

-0.9% 

-2.5% 

-4.3% 

-4.0% 

-2.4% 

Large tower B 

8 .8% 

2.3% 

1.3 % 

-0.5% 

-0.3% 

-4.7% 

-2.4% 

-1.5% 

-6.5% 

- 


Table 3: Relative differences in eight bin energy spectra to evaluate systematic error due to position resolution. Small tower, Large tower A, and 
Large tower B correspond to the rapidity ranges from 10.76 to infinity, from 8.99 to 9.22, and from 8.82 to 8.99, respectively. 


LHCf Energy 
[GeV] 

100 

-500 

500 

-1000 

1000 

-1500 

1500 

-2000 

2000 

-2500 

2500 

-3000 

3000 

-3500 

3500 < 

Small tower 

23.2% 

8.4% 

4.5% 

3.9% 

2.5% 

2.4% 

1.9% 

2 .0% 

Large tower A 

-54.2% 

-1.0% 

0.9% 

-2.5% 

-0.6% 

-3.3% 

-3.6% 

-5.0% 

Large tower B 

-11.0% 

-0.0% 

1.4% 

0 .6% 

-1.1% 

-1.0% 

-3.7% 

-5.7% 





Figure 2: Energy spectra of neutron-like events measured by the Arml and Arm2 detectors. The left panel shows the results from the small towers, 
and the center and right panels show the results for the large towers. The horizontal axes represent the reconstructed energy. The vertical bars 
represent the statistical (they are negligibly small) and systematic uncertainties except the energy scale and the luminosity uncertainties. 
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here because it is treated in the unfolding process dis¬ 
cussed in Section 14.21 In spite of the difference in de¬ 
tector response, data from both Arms show the same 
characteristic feature of the spectra. 

Figures [2 compares the energy spectra measured by 
the Arml detector with the MC predictions. Colored 
lines indicate MC predictions by EPOS 1.99 (magenta), 
QGSJET 11-03 (blue), SYBILL 2.1 (green), DPMJET 
3.04 (red), and PYTHIA 8.145 (yellow). The model 
spectra were obtained from full detector simulations and 
taking account of the same reconstruction process as 
the experimental data. None of the models perfectly 
matches the experimental data. Experimental result in¬ 
dicates the hardest spectra than any other MC predic¬ 
tions. 


4.2. Spectra unfolding 

To estimate the true energy distribution, we used the 
multidimensional-spectra unfolding method Hi with 
the variables energy and transeverse momentum (pr)- 
To create training samples for the unfolding process, we 
used MC simulations with neutrons having a flat energy 
spectrum from 50 to 3500 GeV and an uniform injec¬ 
tion to the detector plane. The training samples were 
reconstructed by using the same method as the experi¬ 
mental data. Performance of the unfolding method was 
checked by applying the unfolding process to the MC 
spectra and comparing the results with the true spectra. 
The upper panel in Figure [4] shows the unfolded spec- 



Figure 4: Comparison of unfolded spectra with true spectra for DPM¬ 
JET 3.04 and EPOS 1.99 models at the small tower of Arml. Bottom 
panel shows the ratio of the unfolded spectra to the true spectra. 


tra for the DPMJET 3.04 and EPOS 1.99 models to¬ 
gether with the true spectra at the small tower of Arml. 
Here ’’true spectra” means the true neutron energy dis¬ 
tribution of the MC events after acceptance and trigger 
threshold were applied. The bottom panel shows the 
ratio of the unfolded spectra to the true spectra. The 
differences between the unfolded and true spectra were 
mostly within 20% except at the highest energy bins 
(50-100%). These systematic differences were due to 
the choice of the flat energy distribution as a training 
sample. We found that the difference did not strongly 
depend among five input models. Thus we applied an¬ 
other correction; dividing the unfolded spectra by the 
average of differences. The differences among the five 
models, typically ±10%, were considered as a part of 
the systematic uncertainties. 

Figure [5] shows the unfolded experimental spectra 
measured by the Arml and Arm2 detectors for each ra¬ 
pidity range. The hatched areas show the Arml system¬ 
atic errors, and the bars represent the Arm2 systematic 
errors. The detection efficiency of neutrons and the cor¬ 
rection in the PID efficiency and purity were also con¬ 
sidered. The results below 500 GeV were not shown 
because of the large systematic errors. The unfolded 
spectra from both Arms show good agreement within 
systematic errors. 

The experimental spectra were combined according 
to the previously used method d. It was assumed 
that the systematic uncertainties caused by the energy 
scale, PID correction, beam center position, multihit 
events, and position resolution have bin-to-bin correla¬ 
tions while the other elements are independent between 
bins. The systematic uncertainties except the unfolding 
processes were thought to be fully uncorrelated between 
Arml and Arm2. Because we treated the systematic un¬ 
certainties of the unfolding processes identical in Arml 
and Arm2, these errors were quadratically added after 
the combining process. 

The differential neutron production cross sections 
dcr n ldE were calculated from the unfolded experimen¬ 
tal spectra by using 


dcr„ldE = 


dN(Ar/, A E) 1 2n 
A E L df’ 


( 2 ) 


where dN(At], AE) means the number of neutrons ob¬ 
served in the each rapidity range. A//, and each energy 
bin, A E. L is the integrated luminosity corresponding to 
the data set. The cross sections are summarized in Ta¬ 
ble [A3] Figure [6] shows the combined Arml and Arm2 
spectra together with the model predictions. The ex¬ 
perimental results indicate that the highest neutron pro¬ 
duction rate compared with the MC models occurs at 
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Figure 3: Measured Arml energy spectra of neutron-like events together with MC predictions. Left panel shows the results for the small tower, 
and the center and right panels show the results for the large tower. The vertical bars represent the statistical (they are very small) and systematic 
uncertainties except the energy scale and luminosity uncertainties. Colored lines indicate MC predictions by EPOS 1.99 (magenta), QGSJET 11-03 
(blue), SYBILL 2.1 (green), DPMJET 3.04 (red), and PYTHIA 8.145 (yellow). 



Energy [GeV] Energy [GeV] Energy [GeV] 


Figure 5: Unfolded energy spectra of the small towers (77 > 10.76) and the large towers (8.99 < 77 < 9.22 and 8.81 < 77 < 8.99). The hatched areas 
show the Arml systematic errors, and the bars represent the Arm2 systematic errors except the luminosity uncertainty.. 





Figure 6 : Comparison of the LHCf results with model predictions at small tower (77 > 10.76) and large towers (8.99 < 77 < 9.22 and 8.81 < 77 < 
8.99). The black markers and gray hatched areas show the combined results of the LHCf Arml and Arm2 detectors and the systematic errors, 
respectively. 
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Table 4: Hadron-to-photon ratio for experiment and MC models. The 
number of neutrons with energy above 100 GeV was divided by the 
number of photons with energy above 100 GeV. The rapidity intervals 
corresponding to the small tower. Large tower A, and Large tower B 
are rj > 10.76, 9.22 > rj > 8.99, and 8.99 > 77 > 8.81, respectively. 


NJNy 

Small 

Large A 

Large B 

Data 

3.05+0.19 

1.26+0.08 

1.10+0.07 

DPMJET 3.04 

1.05 

0.76 

0.74 

EPOS 1.99 

1.80 

0.69 

0.63 

PYTHIA 8.145 

1.27 

0.82 

0.79 

QGSJET 11-03 

2.34 

0.65 

0.56 

SYBILL 2.1 

0.88 

0.57 

0.53 


the most forward rapidity. The QGSJET 11-03 model 
predicts a neutron production rate similar to the experi¬ 
mental results at the largest rapidity range. However, the 
DPMJET 3.04 model predicts neutron production rates 
better at the smaller rapidity ranges. 

The neutron-to-photon ratios (N„/N y ) at three differ¬ 
ent rapidity regions are extracted after unfolding and 
summarized in Table [4] Here N„ and N y are the number 
of neutrons and number of photons, respectively, with 
energies greater than 100 GeV. The numbers of pho¬ 
tons were obtained from the previous analysis 0 and 
the same analysis for the pseudo-rapidity range 8.99- 
9.22 defined in this study. The experimental data show 
the highest neutron ratio compared with the hadronic- 
interaction models for all rapidity regions. 

5. Summary and discussions 

An initial analysis of neutron spectra at the very for¬ 
ward region of the LHC is presented in this paper. The 
data were acquired in May 2010 at the LHC from \[s 
= 7 TeV proton-proton collisions with integrated lumi¬ 
nosity of 0.68 nhr 1 and 0.53 nb _1 for the LHCf Arml 
and Arm2 detectors, respectively. 

The neutron energy spectra were analyzed in the three 
different rapidity regions. The results obtained from the 
two independent calorimeters of Arml and Arm2 are 
consistent each other. The measured spectra were com¬ 
bined and unfolded by using a two-dimensional unfold¬ 
ing method based on Bayesian theory. 

The experimental results, both in folded and unfolded 
were compared with the MC predictions of QGSJET 
11-03, EPOS 1.99, DPMJET 3.04, PYTH1A 8.145, 
and SYBILL 2.1; however, no model perfectly repro¬ 
duces the experimental results in the whole pseudo¬ 
rapidity range. Moreover, compared with the hadronic- 
interaction models, the experimental results show the 


most abundant neutron production relative to the pho¬ 
ton production rate. 

Because of the difference in the phase space, neutrons 
at lower rapidity carry more energy than those at higher 
rapidity. EPOS and QGSJET 110, current standard mod¬ 
els in the air shower analysis, underestimate the neutron 
production to about 30% at the lower rapidity regions. It 
is interesting to investigate how these differences affect 
to the development of air showers. 

The differential cross sections for neutron production 
at very forward rapidity were measured at the Intersect¬ 
ing Storage Ring for proton-proton collisions at \[s — 
30.6-62.7 GeV j 190 and by the PHENIX experiment 
at the Relativistic Heavy Ion Collider for proton-proton 
collisions at \[s — 200 GeV jTJ] . The results of previ¬ 
ous experiments were consistent with Feynman x (xf ) 
scaling and do not depend on the collision energy. To 
accurately extrapolate the hadronic-interaction models 
to the high-energy range (independently of whether or 
not the xf scaling in the neutron production is still rel¬ 
evant at energies in the TeV range) is also an important 
issue. LHCf will extend the energy range to test the 
Feynman scaling with the proton-proton collision data 
obtained at yfs- 0.9, 2.76, 7 TeV and to be obtained at 
V^=13TeV in 2015. 
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Table A.5: Differential neutron production rate dcr„/dE [mb/GeV] for each rapidity range. 
Cross section [mb/GeV] 


Energy [GeV] 

Small tower (77 > 10.76) 

Large tower A (8.99 < 77 < 9.22) 

Large tower B (8.81 < 77 < 8.99) 

500-700 

(s^i+o.s^xio - 3 

(5.91±0.66)xl0~ 4 

(6.09±0.84)xl0~ 4 

700-900 

(7.48+0.95)xl0~ 5 

(6.38±0.82)xl0~ 4 

(7.15±0.67)xl0~ 4 

900-1100 

(9.32±1.54)xl0~ 5 

(7.70±0.66)xl0~ 4 

(8.60±0.92)xl0~ 4 

1100-1300 

(1.26±0.19)xl0~ 4 

(8.41±0.75)xl0~ 4 

(9.42±1.05)xl0~ 4 

1300-1500 

(1.58+0.32)xl0- 4 

(9.18±0.80)xl0~ 4 

(9.71±0.91)xl0~ 4 

1500-1700 

(1.95±0.31)xl0- 4 

(9.54±0.92)xl0- 4 

(9.34±0.59)xl0 - 4 

1700-1900 

(2.45+0.35)xl0 ^ 4 

(9.03±0.68)xl0- 4 

(8.12±0.78)xl0- 4 

1900-2100 

(2.72±0.36)xl0- 4 

(8.21±0.81)xl0- 4 

(7.19±0.68)xl0- 4 

2100-2300 

(2.92±0.27)xl0- 4 

(6.90±0.82)xl0- 4 

(5.81±0.55)xl0- 4 

2300-2500 

(2.98±0.28)xl0- 4 

(6.17±0.52)xl0- 4 

(4.25±0.53)xl0~ 4 

2500-2750 

(2.98±0.34)xl0~ 4 

(4.21±0.44)xl0 ~ 4 

(2.94±0.54)xl0 ~ 4 

2750-3000 

(2.82±0.48)xl0- 4 

( 2 . 20 ± 0 . 68 )xl 0~ 4 

(1.39±0.65)xl0~ 4 

3000-3250 

(2.49+0.78)xl0~ 4 

(1.10±0.55)xl0 ~ 4 

(6.07±3.39)xl0- 5 

3250-3500 

(2.32±1.06)xl0~ 4 

(2.45±1.70)xl0~ 5 

(5.75±3.76)xl0~ 6 
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